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We show that the interplay between geometric criticality and dynamical fluctuations leads to a 
novel universality class of the contact process on a randomly diluted lattice. The nonequilibrium 
phase transition across the percolation threshold of the lattice is characterized by unconventional 
activated (exponential) dynamical scaling and strong GrifBths effects. We calculate the critical 
behavior in two and three space dimensions, and we also relate our results to the recently found 
infinite-randomness fixed point in the disordered one-dimensional contact process. 
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Nonequilibrium systems can undergo continuous phase 
transitions between different steady states. These tran- 
sitions are characterized by collective fluctuations over 
large distances and long times similar to the behavior of 
equilibrium critical points. Examples can be found in 
population dynamics and epidemics, chemical reactions, 
growing surfaces, and in granular flow and traffic jams 
(for recent reviews see, e.g., Refs. 

If a nonequilibrium process is defined on a randomly 
diluted spatial lattice, its dynamical fluctuations coex- 
ist with geometric fluctuations. Site or bond dilution 
defines a percolation problem for the lattice with a geo- 
metric phase transition at the percolation threshold 0. 
In this Letter we address the question of how the inter- 
play between geometric criticality due to percolation and 
dynamical fluctuations of the nonequilibrium process in- 
fluences the properties of the phase transition. 

Our starting point is the contact process |^, a proto- 
typical system exhibiting a nonequilibrium phase transi- 
tion. It is defined on a d-dimensional hypercubic lattice 
{d > 2). Each site can be active (occupied by a particle) 
or inactive (empty). In the course of the time evolution, 
active sites infect their neighbors, or they spontaneously 
become inactive. Specifically, the dynamics is given by 
a continuous-time Markov process during which particles 
are created at empty sites at a rate Xn/(2d) where n is 
the number of active nearest neighbor sites. Particles are 
annihilated at unit rate. For small birth rate A, annihila- 
tion dominates, and the absorbing state without any par- 
ticles is the only steady state (inactive phase). For large 
birth rate A, there is a steady state with finite particle 
density (active phase). The two phases are separated by 
a nonequilibrium phase transition in the directed perco- 
lation 0,01 universality class at some A = A°. 

We introduce quenched site dilution by randomly 
removing lattice sites with probability p. The result- 
ing phase diagram of the site-diluted contact process is 
sketched in Fig. ^ For small impurity concentrations 
below the percolation threshold of the lattice, p < Pc, 
the active phase survives, but the critical birth rate in- 
creases with p (to compensate for the missing neighbors) . 
Right at the percolation threshold the active phase sur- 



vives on the infinite percolation cluster for A > A*. The 
(multi) critical birthrate A* must be smaller than the crit- 
ical birthrate of the one-dimensional (ID) contact process 
because the critical percolation cluster is connected, in- 
finitely extended, and its fractal dimension is Df > 1. 
For p > Pc, no active phase can exist because the lattice 
consists of disconnected clusters of finite size that do not 
support a steady state density of active sites. 

The contact process on a site-diluted lattice therefore 
has two nonequilibrium phase transitions, separated by a 
multicritical point. For p < pc, the transition (marked by 
"a" in Fig.nj is expected to be in the universality class of 
the generic disordered contact process 0, 0, IT3 w hich 
has reattracted considerable attention recently '15, 
In contrast, the phase transition across the percolation 
threshold Pc of the lattice for A > A* (transition "b" in 
Fig. ^ has received much less attention. 

In this Letter, we show that the interplay between 
geometric criticality and dynamic fluctuations leads to 
a novel universality class for this nonequilibrium phase 
transition. Even though the transition is driven entirely 
by the geometry of the lattice, the dynamical fluctua- 
tions of the contact process enhance the singularities in 
all quantities involving dynamic correlations. Our results 
can be summarized as follows. The dynamical scaling is 
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FIG. 1: Schematic phase diagram of a site diluted contact 
process as function of impurity concentration p and birth rate 
A. There is a multicritical point at p = pc and A = A*. The 
phase transition (b) across the percolation threshold of the 
lattice is the topic of this Letter. 
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not of conventional power-law form but activated, i.e., 
the relation between correlation length and correla- 
tion time is exponential. 



(1) 



with the critical exponent ip being equal to the fractal 
dimension of the critical percolation cluster, ip — Df. As 
a result, the long-time decay of the density p of active 
sites at p = pc is ultra-slow. 
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The exponent S = PdivcDj) is determined by Df to- 
gether with the order parameter and correlation length 
exponents, Pc and i^c, of the lattice percolation transi- 
tion In contrast to the enhanced dynamical singu- 
larities, the exponents of static quantities like the steady 
state density pst and the spatial correlation length are 
identical the corresponding lattice percolation exponents. 
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where A — p — pc measures the distance from the per- 
colation threshold. Off criticality, i.e., away from the 
percolation threshold, we find strong Griffiths effec ts jlE 
characterized by a non-exponential density decay [l^ 
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where the nonuniversal exponents z' and z" diverge as 
p ^ Pc- In the remainder of the Letter 
we sketch the derivation of these results and calculate ex- 
ponent values and additional observables. We also relate 
our results to the disordered ID contact process 0, 
and to the diluted quantum Ising model jl9|. 

Let us start by considering the steady state density 
Pst of active sites (i.e., the order parameter of the transi- 
tion). A nonzero steady state density can only develop on 
the infinite percolation cluster; finite clusters do not con- 
tribute because they eventually go into the inactive state 
via a rare fluctuation. For A > Ac, the infinite cluster 
is in the active phase. The total steady state density is 
proportional to the number of sites in the infinite cluster, 
Pst ~ Pooiv) ~ iPc — p)^" ■ The order parameter expo- 
nent of the nonequilibrium transition is therefore identi- 
cal to that of the lattice percolation transition, (3 = f3c, as 
stated in 0. To determine the spatial correlation length 
we note that the correlations of the contact process 
cannot extend beyond the connectedness length of the 
percolating lattice because sites on different percolation 
clusters are decoupled. On the other hand, for A > A*, all 
sites on the same cluster are strongly correlated in space 
even though they collectively fluctuate in time. We thus 
conclude ^ ~ and = i^c in agreement with Q . 



We now study the time dependence of the density p{t) 
of active sites, starting from a completely active lattice. 
We first consider the contact process on a single per- 
colation cluster of finite size (number of sites) s. For 
A > A* such a cluster is locally in the active phase. It 
therefore has a metastable state with a nonzero density 
of active sites. This metastable state can decay into the 
inactive state only via a rare collective fluctuation involv- 
ing all sites of the cluster. The probability for such a rare 
event decreases exponentially with the size s of the clus- 
ter. Therefore, the life time ts of the metastable active 
state on a cluster increases exponentially with its size s. 
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where to is some microscopic time scale. The prefactor 
A{\) vanishes at the multicritical point, A{\^,) = 0, and 
increases with increasing A. The number of sites s of a 
percolation cluster is connected to its linear size Rs via 
s ~ i?f ^ . Therefore, establishes the exponential rela- 
tion between length and time scales, Intg ^ Rs leading 
to activated dynamical scaling 

After having analyzed a single cluster we turn to the 
full percolation problem. From classical percolation the- 
ory, we know that close to Pc the number of occupied 
clusters of size s per lattice site (excluding the infinite 
cluster for p < pc) obeys the scaling form 
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The scaling function f{x) behaves as 
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where i3i,i?2 are constants. The exponents Tc and Cc 
determine all critical exponents of the percolation tran- 
sition of the lattice including the correlation length ex- 
ponent Vc — (tc — ^)/id'0'c), the order parameter ex- 
ponent f3c = (tc — 2)/ac, and the fractal dimension 
Df — d/ {tc ~ 1) of the percolating cluster (jj. 

In order to obtain the total density of active sites for 
the contact process on the diluted lattice, we sum the 
number of active sites over all percolation clusters. Com- 
bining the cluster size distribution JH)) with the lifetime 
of the metastable active state (O leads to 



pit, A) ~ J ds s nsiA)exp[-t/tsis)] (12) 
Right at the percolation threshold, this reduces to 

p{t,0) ^ I ds si-^'^exp[-(t/toe'^")] . (13) 



The leading behavior of this integral can be found by 
noticing that only islands with size s > s„iin(i) = 
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A ^ ln(t/to) contribute at time t. The critical long-time 
dependence of the total density is thus given by 
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This completes the derivation of 10) with the critical ex- 
ponent S given by (5 = Tc — 2 = Pr J (vrDf) in agreement 
with general scaling arguments [l^ |la| ■ 

We now consider the behavior of the density off criti- 
cality. In the inactive phase, p > Pc, the time dependence 
of the density is given by 

p(t,A)- y dssi-^-exphSisAi/'^- - (i/ioe^")]. (15) 

For long times, the leading behavior of the integral can 
be calculated using the saddle-point method, giving 
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equivalent to The nonuniversal exponent z' is given 
by z' = (Ad/Bi)A-i/-c 

In the active phase, p < Pc, there is a nonzero steady 
state density p^t coming from the infinite percolation 
cluster. However, the approach of the density towards 
this value is still determined by the slow decay of the 
metastable states of the finite percolation clusters 
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Using the saddle point method to calculate the leading 
long-time behavior gives (for p < pc) 
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This completes the derivation of © with z" = 
{Ad/ 32)1^1-^/"-= ^ The non-exponential off- 

critical relaxation of the density (|5I6(I is characteristic of 
a Griffiths region in the contact process 0, . We also 
point out that time and spatial correlation length enter 
these equations in the form of the combination ln(t)/^^^ 
again characteristic of activated scaling. 

We now turn to the influence of an external source 
field h that describes spontaneous particle creation at a 
rate h at each lattice site. To determine the steady state 
density as a function of h we again start by considering 
a single percolation cluster of size s. For A > A*, the 
cluster is active if at least one particle has been spon- 
taneously created on one of the s sites within the life 
time ts{s) = toe"^^ ■ For small h, the average number 
of particles created on a cluster of size s within time tg 
is Ms{h) = hsts — hstoe^'^ . If Alg > 1, the cluster is 
(almost) always active. If Ms < 1, it is active with a 
probability proportional to Ms- The total steady state 
density is obtained by summing over all clusters 
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5/36 


0.417 
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0.875 


= Df = d — Pc/vc 


91/48 


2.523 


5 = l3a/{VcDf) 


5/91 


0.188 



TABLE I: Critical exponents of the nonequilibrium phase 
transition at p = pc in two and three space dimensions. 



Evaluating this integral analogously to the time- 
dependent density H12|l yields, for small fields h, 
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where ho l/^o- P — Pc, the relation between density 
and field is logarithmic, as expected from activated scal- 
ing. Off criticality, we find strong Griffiths effects similar 
to those in the time-dependence of the density. 

The above results can also be derived from a scaling 
theory. In the active phase, the density is proportional 
to the number of sites in the infinite percolation cluster. 
Thus, its scale dimension must be Pd^c- Time must 
enter via the scaling combination h\{t)h^^ reflecting the 
exponential dependence of the life time (UJ on the cluster 
size. The field h, being a rate, scales like inverse time. 
We therefore obtain the following scaling form: 

p[A,ln(t),ln(l//i)] = b'^-'/"" plAb-^'"-' ,\n{t)b'l' ,\iv{l/h)V''] 

(23) 

where b is an arbitrary (length) scale factor and tp = D f. 
This form is consistent with all our explicit results. 

All critical exponents of the nonequilibrium phase 
transition are determined by the classical percolation ex- 
ponents of the lattice. In two space dimensions, their 
values are known exactly and in three dimensions they 
are known numerically with high accuracy 0. Table I 
shows numerical exponent values for these cases. 

We also briefly comment on the early time behavior. 
For A > A* each percolation cluster is locally in the active 
phase. Starting from a single active seed, the cloud of ac- 
tive sites thus initially grows ballistically, i.e., the radius 
of the cloud grows linearly with time, until a metastable 
state is reached in which it covers the entire percolation 
cluster. The time required for this initial spreading on an 
island of size s is ti{s) ^ Rs ^ s^^^f . As discussed above, 
the metastable state decays only at the much larger time 
scale ts{s) ~ e^'*. We thus arrive at the somewhat sur- 
prising conclusion that the early time behavior of the 
contact process on our diluted lattice is much faster than 
the logarithmically slow long-time decay of the density. 

In the remaining paragraphs, we discuss the generality 
of our results, compare them to the transition in the di- 
luted quantum Ising model |l9l | and to the recently found 
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infinite-randomness critical point in a random ID contact 
process 0, . We also compare to a general classifica- 
tion of phase transitions with quenched disorder . 

The logarithmic time and field dependencies l(2Jl and 
(|20|l at the nonequilibrium phase transition at p — Pc 
as well as the strong Griffiths effects in its vicinity are 
the direct result of combining the spectrum of percola- 
tion cluster sizes (jHJ with the exponential dependence 
(0) of the life time on the cluster size. We therefore 
expect similar behavior in other diluted equilibrium or 
nonequilibrium systems that share this exponential rela- 
tion between length and time scales. One example is the 
diluted transverse field Ising model In this system, 
the energy gap of a cluster decreases exponentially with 
its size. As a result, the scaling behavior at the quantum 
phase transition across the percolation threshold of the 
lattice is very similar to the one found in this paper. 

Recently, the critical point of the ID contact pro- 
cess with spatial disorder was found to be of infinite- 
randomness type 0, 0|. Analogous behavior is ex- 
pected for the generic disordered directed percolation 
transition in higher dimensions, e.g., the transition at p < 
Pc in the site-diluted contact process (transition a in Fig. 
. Our critical point at the percolation threshold shares 
some characteristics with these infinite-randomness criti- 
cal points, notably the exponential relation between cor- 
relation length and time as well as the logarithmically 
slow decay of the total density. However, it belongs to a 
different universality class with novel critical exponents. 
Moreover, the early time behavior is different (logarith- 
mically slow at the generic infinite randomness critical 
point but of power- law type at our transition) . 

Lastly, we point out that our results are in agreement 
with a general classification of phase transitions with 
quenched disorder (and short-range interactions) accord- 
ing to the effective dimensionality dg// of the droplets 
or clusters [i^. Three cases can be distinguished: (i) If 
the clusters are below the lower critical dimension of the 
problem, de// < d'^ ^ the critical behavior is of conven- 
tional power-law type and the Griffiths effects are expo- 
nentially weak, (ii) If def j = d~ , the critical point shows 
activated scaling accompanied by strong, power-law Grif- 
fiths effects. This case is realized in random transverse 
field Ising magnets 0, |^ as well as in our diluted con- 
tact process, (iii) If deff > d^ , the phase transition is 
smeared because locally ordered clusters can undergo the 
phase transition independently from the bulk. This oc- 
curs, e.g., for some metallic quantum magnets J22| or for 
the contact process with extended defects UM- 

In conclusion, we have shown that the contact process 
on a diluted lattice has unusual properties. The inter- 
play between geometric criticality and dynamical fluctu- 
ations leads to a novel universality class with activated 
scaling and ultraslow dynamics. Interestingly, despite its 
ubiquity in theory, experimental observations of directed 
percolation scaling j2^] are very rare. Our results suggest 



that peculiar disorder effects may be responsible for this 
in at least some of the experiments. 
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